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Abstract 

A new potential energy surface for the electronic ground state of the simplest triatomic anion 
is determined for a large number of geometries. Its accuracy is improved at short and large 
distances compared to previous studies. The permanent dipole moment surface of the state is 
also computed for the first time. Nine vibrational levels of H3 and fourteen levels of D3 are 
obtained, bound by at most ~ 70 cm -1 and ~ 126 cm -1 respectively. These results should guide 
the spectroscopic search of the Hg" ion in cold gases (below 100K) of molecular hydrogen in the 
presence of ions. 

PACS numbers: 82.20.Kh, 34.20.-b, 82.30.Fi 
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I. INTRODUCTION 



Collisions involving hydrogen atoms, molecules, and their positive (H + , H 2 ) and negative 
(H~) ions play an important role in chemistry and evolution of neutral or negatively-charged 
hydrogen plasma such as laboratory hydrogen plasma, the interstellar medium (ISM), at- 
mospheres of the Sun and other stars [1] as well in the Earth atmosphere. Binary collisions 
of these species are simple enough to be treated using first principle methods without any 
adjustable parameter. Therefore, such processes are often used as benchmarks for testing 
theoretical methods. 

The existence of bound states of the ion in a linear configuration for the three nuclei 
has been suggested in 1937 j^j, but has been questioned since then. Different ab initio 
calculations pHl 1 1 have been giving contradictory results for the depth of the van der Waals 
well. The estimated error bar in the calculations was comparable or larger than obtained 
binding energies, so the calculations could not predict for sure if the Hg ion is stable. 
Although first observations of ions have been reported in low resolution experiments as 



early as in 1974 



10 



12 



17] , they could not warrant either the stability of H 3 . It is only in 
that theory became precise enough to confirm the stability of Hj bound 



the 1990's 
states. 

At present, there are three available accurate potential energy surfaces (PES) of H^": 
The PES by Starck and Meyer [loj], by Belyaev and Tiukanov 19], and by Panda and 
Sathyamurthya [llj]. In the following, these three PES will be referred to as PES-SM, 
PES-BT, and PES-PS respectively. In addition, Belyaev, Tiukanov and others 19H24] have 
obtained PES of excited electronic states of Hg and their non-Born-Oppenheimer couplings 
with the ground state. The excited electronic states are unstable with respect to electron 
autodetachment. Bound state calculation based on the PES-SM have been performed for 



H7, H 2 D~, and D 2 H~ in Ref. 



lOj and using the PES-BT for H 3 in Ref. [l6(. It is worth to 



mention a stand alone study by Robicheaux 18] confirming the stability of (and similar 
anions), where the PES was derived from the scattering length for the electron-H 2 collisions, 
and from the polarizability of H 2 . 

While never observed in the ISM, the collisions between H 2 and H~, and of their isopoto- 



logues, have been studied in a number of laboratory experiments back to the 1950 's 
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Such collisions should also play a role in processes in tokamaks, especially if the negative 
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ion source (D~) is used for the tokamak neutral beam injectors 32|, |33j. Many studies have 
been devoted to the calculation of elastic and inelastic cross-sections for H 2 -H~ collisions 
and for all isotopologues, [ll, [2^, 34 - 45 1 , most of them for energies significantly larger than 



1 eV above the lowest dissociation limit H 2 + H~. However, at such energies the theoret- 
ical cross-sections may not be reliable: the dissociation energy of the H 3 ground state lies 
only 0.75 eV above the H3 dissociation energy, and the non-Born-Oppenheimer interaction 
between the ground states of H3 and could be significant. The non-Born-Oppenheimer 
interactions have been taken into account only in Ref. 46] within a reduced-dimensionality 
approach, where only linear geometries of have been taken into account. Finally, the 
photodissociation of H 2 D~ has been studied in a simplified approach based on the Franck- 
Condon overlap between a bound level and a scattering state of H 2 D~ (4^], since the ab 
initio dipole moments of H3 were not available at that time. 

The present study is mainly motivated by the formation of bound states in low-energy 
collisions between H 2 and H~. In such collisions Hj can be formed only if a third body (other 
than H 2 or H - ) participates (three-body recombination - TBR) or if a photon is emitted 
(radiative association - RA). Both processes could be relevant for the chemistry of cold 
interstellar clouds, if H~ is present |48] . Note that the H~ ion has not been detected so far 
in the interstellar medium (ISM): it cannot be directly observed by usual photoabsorption 
spectroscopy because H~ has only one bound electronic state. 

The evaluation of the cross-section for radiative association of H 2 and H at low energy 
(10-30K) requires an accurate PES with a precision around 1 cm -1 or better, and permanent 
dipole moment surfaces (PDMS) for Hj. In the present study we calculated a new ab initio 
PES for the H3 electronic ground state on a dense and large grid for internal coordinates, 
using a larger electronic basis set than those of PES-SM and PES-PS calculations. Special 
care is taken in order to account for the long-range behavior of the surface. We also obtain 
for the first time the PDMS of the electronic ground state. We constructed Fortran 
subroutines that calculate PES and PDMS values for any arbitrary geometry using B-spline 
interpolation procedures, and we determine the bound states of Hj and D^. The RA reaction 
will be treated in a separate study. 

The article is organized in the following way. In the next section, we discuss the calcu- 
lation of the new PES and the interpolation procedure. In Section [TTT1 we compare the new 
PES with the PES from previous studies. Section [TV] is devoted to the calculation of bound 
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levels of Hj and and Section IVl presents our results on the dipole moment of Hj. Section 
IVII is the conclusion. Atomic units (a.u.) for distances (1 a.u. = 0.0529177 nm) and for 
energies (1 a.u. = 219474.63137 cm -1 ) will be used throughout the paper, except otherwise 
stated. 



II. AB INITIO CALCULATION AND INTERPOLATION OF THE GROUND 
STATE POTENTIAL SURFACE 



As in Ref. 



10J, we used the coupled-electron pair approximation (CEPA-2) method 49] , 



which is part of the Molpro package 50] . It is a non- variational variant of the configuration 



interaction method for closed-shell molecules. Here we used a considerably larger basis set, 
AV5Z with spdfg basis functions from the Molpro basis library, and a much larger number 
of geometries, than in Refs Q, [ll] • As Hg is a van der Waals molecule, we defined a 3D-grid 
in Jacobi coordinates: r- the distance between two protons, R- the distance from the center 
of mass of the two protons to the third proton, and 7 - the angle between vectors R and r. 
The grid in r is uniform from r = 0.8 a.u. to r = 2.4 a.u. and changes by a step Ar = 0.2 
a.u. The grid in 7 changes from to 90° by a constant step of A7 =15°. The grid in R was 
chosen denser for small R than for large R: the grid points Ri were calculated according to 
Ri = 1.5 + 0.452 exp(i/10) (in a.u.) with i = 1,2, •• • ,48, which makes R changing from 
1.9995 a.u. to 56.4227 a.u. Therefore, the calculations were performed for 9 x 7 x 48 = 3024 
geometries. Notice that the above grid starts at a smaller value of r = 0.8 a.u. than in 



Refs. 



10 



11] (r = 1 a.u.). Indeed, we found that a grid starting at r = 1 cannot properly 
represent the repulsive part of the H 2 ground state potential curve (at fixed R and 7) and 
gives an appreciable error in energies obtained of the lowest H2 vibrational levels and, as a 
result, a comparable error in dissociation energies of — kH 2 (f , j)+H~ when R — > 00. We 
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also used a denser grid in 7 and a longer and denser grid in R than in Refs. 

For each geometry we calculated the potential energy of the ground state, and 
the components of the PDM vector with respect to the principal axes of inertia. The 
obtained ab initio PES and PDMS were used to prepare Fortran subroutines that cal- 
culate PES and PDMS for any arbitrary geometry. Inside the three-dimensional box 
r G [0.8; 2.4], R G [1.9995; 56.4227], 7 G [0; 360°] (central region in Fig. d)) the proce- 
dures interpolate the surfaces using the three-dimensional B-spline method. Outside of the 
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FIG. 1: The five regions of the configuration space in Jacobi coordinates R, r defined for the PES 
calculations, for all values of 7 (see the text). Notice that R is given on a logarithmic scale. 
The central rectangle represents the region of calculated ab initio points, where a 3D B-spline 
interpolation is used. Different methods were used for the extrapolation of the PES and PDMS 
in regions I- IV (see text). The arcs of a circle (dashed lines) represent lines of fixed hyperradii p 
(see Section |IV]) for two values, p% = 2 a.u. and p2 = 8 a.u. The X, Y, and Z axes (needed to 
specify the components of the PDM vector) are the axes of principal moments of inertia. Their 
orientation is indicated in the figure. The Y-axis is out the plane of Hg\ 

box, the procedures rely on analytical formulas for the extrapolation of PES and PDMS. 
To simplify the description of this region, we divided the r x R configuration space in four 
different parts (I- IV) surrounding the ab initio rectangle (Fig. [T]). Regions 11 ab initio calcu- 
lations'''' and II are the only ones relevant to bound and scattering states of the system with 
energies < 2eV above the fi2(f = 0, j = 0)+H~ dissociation. Therefore, the extrapolation 
procedure for the PES and PDMS in the two regions should be physically justified. In con- 
trast, wave functions of bound and continuum states with such energies vanish in regions I, 
III, and IV, and empirical formulas will be used. As it will be discussed in Section IIV[ the 
extrapolation is needed to map the configuration space in Jacobi coordinates on the space of 
hyperspherical coordinates, which will be used for bound state and scattering calculations. 
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In region II, we represent the long-range (in R, at fixed r and 7) potential Vlr for the 
interaction between H2 and H~ as [511 ] : 

V LR (R;r,j) = D as (r) + -JL + -± with (1) 
Cf = -Q(r)P 2 (cos 7 ) and Cf = - [or (r) + a 2 (r)P 2 (cos 7 )] /2 , 

where the first term D as (r) is the sum of H~ and H 2 (r) energies at a given internuclear 
distance r of the H 2 molecule. The second term is the interaction between the electric 



charge of n with the quadrupole moment Q(r) of H 2 (taken from Ref. |52[]), and the third 
term is the interaction of the dipole moment of H 2 induced by H~, involving the second 
order Legendre polynomial P 2 (cos 7 ). The functions ao(r) and a 2 (r) are the isotropic and 



anisotropic polarizabilities of H 2 , for which we used the analytical functions given in Ref. 
that were obtained by fitting the numerical values from Ref. 



531 ] . The dispersion energy 



varying as 1/R 6 and other smaller terms are neglected, which is a good approximation 
because the long-range expansion is only used for R > 56.4 a.u. 

In region I we used the following extrapolation formula in r for fixed R and 7: 

V SR (r ] R, 1 ) = a(R, 1 )e- b ^ r , (2) 

where a(R, 7) and b(R, 7) are functions of R and 7 that are obtained considering the two 
ab initio energies V{r = 0.8; R, 7) and V(r = 1;R, 7) calculated at first two values of the 
coordinate r. In this way, we obtain the quantities a and b given on a two-dimensional grid 
of points in the (P, 7 ) space. Then we used the 2D B-spline interpolation to obtain smooth 
two-dimensional functions a(R, 7) and 6(P, 7). 

In region III, we extrapolate the PES in r and at fixed R and 7 using a dispersion- like 
expression: 

V LR (r; R, 7) = D (R, 7) - , (3) 

where the D {R^) and C 6 (R, 7) (always positive) coefficients are obtained in a way similar 
to the coefficients a and b, considering the two last points V(r = 2.2 a.u.; R, 7) and V(r = 2.4 
a.u.; R, 7). They are also interpolated using the 2D B-spline method for arbitrary values of 
R and 7. 

The behavior of the PES in region IV is described by the short-range (in R) repulsive 
expression at given values of r and 7: 

V SR (R;r, 1 ) = A(r, 1 )e- B ^ R , (4) 



The A(r, 7) and B(r,j) coefficients are obtained in a way similar to the one described 
above, considering the first two values of the ab initio potential energy V(R = 2 a.u.;r, 7) 
and V(R = 2.052 a.u. ; r, 7), and are further interpolated for any (r, 7) using the same 2D 
B-spline procedure. 




r\ (a.u.) r\ (a.u.) 

FIG. 2: Energy contour plots of the electronic ground state PES in the space of the two 
internuclear distances (r\ and r 2 ) for two values of the bonding angle 5 defined in each plot. 
Successive contours differ by 0.006 a.u. The dissociation energy D as = —1.7016828 a.u. of the 
PES corresponds to the infinite separation between H2(r e =1.403 a.u.) and H - (see Table H}. 
The energy minimum (indicated by M) is obtained for a linear configuration, and the top of the 
potential barrier is labeled with B. The lowest contour corresponds to the energy of —1.7034 a.u. 
(left plot) and -1.6974 a.u. (right plot). 

We show in Fig. [2] the H3 PES as a function of two internuclear distances, T\ and r 2 
for two values (180 and 90°) of the bonding angle 5. These coordinates are convenient to 
describe main features of the PES. As expected, the figure is symmetric with respect to the 
7*1 <H- r 2 exchange, and shows the van der Waals potential well in each coordinate (with the 
minimum labeled with M), separated by a potential barrier (with the maximum labeled 
with B) for the exchange of two identical nuclei. As in previous studies, lowest energy is 
found for a linear configuration. The position and energy of M and B are reported in Table 
HI together with other characteristic parameters of the PES. 

The Fortran subroutines calculating the PES and PDMS from the present ab initio data 
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Quantity 


Present study 


Ref. [10] 


Ref. [11] 


Ref. [19] 


D%\ (a.u.) 


-1.701683 


-1.70095 


n/a 


n/a 


r^(a.u.) 


1.403 


1.40 


n/a 


n/a 


£ c ) (a.u.) 


-1.68509 


-1.68562 


n/a 


n/a d ) 


position of B 
ri, r 2 (a.u.), 5 = 
r, i?(a.u.), 7 = 


1.996, 1.996 
1.996, 2.994 


1.997 , 1.997 
1.997, 2.996 


1.999, 1.999 
1.999, 2.999 


1.74, 1.74 
1.74, 2.61 


M £ ) (a.u.) 


-1.703511 


-1.70270 


n/a 


n/a 


position of M 
r, R (a.u.), 7 = 


1.421, 6.069 


1.416, 6.183 


1.419, 5.915 


n/a 


B — M (cm" 1 ) 


4042.9 


3748.62 


3786.42 


n/a 


D as — M (cm" 1 ) 


401.2 


384.05 


384.27 


443.60 


Energies obtained in calculation made separately for H2 and H 


E R 2 ( r =lA) 


-1.174252 


-1.17368 


n/a 


n/a 


En- 


-0.527429 


-0.52727 


n/a 


n/a 



TABLE I: Different quantities characterizing the Ho ground state PES obtained in the present 



study compared with previous calculations 



id, 
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19( | . a) asymptotic energy at infinite separation 



between H2(r e ) and H ; b) internuclear distance corresponding to the minimum of the H2 dimer 



191 ] gives the height of 



potential; c) energy of the maximum of the barrier, B in Fig. [2 d) Ref. 
the barrier (0.624 eV) with respect to the P±2(v = 0,j = 0)+H _ dissociation; e) energy of the PES 
minimum, M in Fig. [2j 



using the interpolation and extrapolation procedures described above are available with the 
journal's EPAPS service and can be provided by the authors upon a request. 



III. COMPARISON WITH PREVIOUS STUDIES 

Figures El HI and [5] illustrate the comparison of the PES in Jacobi coordinates obtained 
in the present study with the results of Refs. Q, Q for three values of 7 =0°, 30°, and 90° 
respectively. Each figure gives the PES for one value of 7 and eight values of r versus the 
coordinate R. The origin of potential energy yielded by the ab initio procedure corresponds 



S 



to an infinite separation of all electrons and nuclei. Because the absolute energy of the 



PES of Ref. 



llj is unknown, the origin of this surface is chosen in such a way that the 



asymptotic energy of the infinite separation between H 2 (r =1.4 a.u.) and H is the same as 
in the present study. 

-1.63 
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FIG. 3: (Color online) Comparison of the present H 3 PES (circles) in Jacobi coordinates with 



n 



111 ] (solid blue lines with no 



the results from Ref. [10( (crosses on the dashed lines) and Ref. 
symbols) for 7 = 0°. The PES are shown as a function of R for eight different values of r. The 
r = 0.8 a.u. curve is not shown here because it is located at a higher energy. Symbols indicate the 
geometries for which the actual ab initio calculations have been performed. The continuous lines 
connecting the points are obtained by interpolation. No ab initio energies are available for the PES 



of Ref. 



Ill ]: the solid lines without symbols are obtained directly from the analytic formula and 



n 



the Fortran subroutine provided in Ref. ll|. Notice the logarithmic scale in R variable. The inset 
enlarges the region of the minimum of the PES close to equilibrium region at r = 1.4 a.u. 
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FIG. 4: (Color online) Same as Fig. [3] except 7 = 30°. 
The insets of these figures, together with the data in Table H demonstrate that the present 

n 

ab initio energies are systematically lower that those of Ref. [10|] for all geometries. For 
instance, the potential well depth has increased by 0.01% (around 160 cm -1 ) at r e = 1.4 



a.u. This was expected, as we use the same CEPA-2 method than in Ref. [10] with a larger 
basis set, leading then to a PES with an improved accuracy. The overall behavior of the 
present PES and PES-SM as a function of R is similar at small and large distances. In 
contrast, a noticeable difference of PES-PS with the present results and the PES-SM is the 
presence of a potential well at 7 = 90° near R = 7.6 a.u. (see the inset in Fig. [5]) with a 
significant depth (~90 cm -1 for r = 1.4 a.u.). For 7 = 90° the C3/R 3 long-range energy 
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FIG. 5: (Color online) Same as Fig. [3] except 7 = 90°. The two insets, (b) and (c), show the 
behavior of the potentials for r = 1.4 a.u. on two different magnification scales: Region (b) is the 
magnified version of (a), region (c) is magnified (b). 

contribution is positive and C4/-R 4 is negative (see Eq. []]). Therefore, their sum combined 
with the short-range interaction produces a potential curve (for fixed r and 7 near 90° ) that 
has a pronounced minimum in PES-PS and just a should-like feature in the present PES 
(see inset (c) in Fig. [5]). The curve is repulsive asymptotically. A possible reason for the 
above differences is the somewhat smaller basis set (d-aug-cc-PVTZ) used in calculation of 



PES-PS 



The numerical data available about PES-BT 19j,|23j is given in Tabled! Since the calcula- 
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tion method (diatomics-in-molecule -DIM) explicitly uses the ab initio data (wave functions 
and energies) obtained for fragments (H 2 , H~, , H 3 ) in separate ab initio calculations, it 
gives automatically the correct dissociation energy. 

Since we have calculated PES at relatively large distances R, we can extract ab initio 
values (labelled with "ai") of the corresponding long-range coefficients Cf 1 and C4 1 and 
compare them with the theoretical values Cf 1 and Cf 1 of Eq. ([1]). We first make a fit of 
C^ 1 to the ab initio data by fixing Cf 1 to the theoretical value C| . In this way the value of 
C4 1 agrees with the value of Cf 1 within 0.3% for 7 = 0, and 5 % for 7 = 90°at r = 1.4 a.u. 
Then, we did the opposite: we fitted Cf 1 to the ab initio data by fixing Cf 1 to Cf 1 . In this 
case we found that fitted and theoretical values of C3 agree within 0.1 % for 7 = and 2 % 
for 7 = 90°at r = 1.4 a.u. The PES-PS differs significantly from the analytical long-range 
behavior of Eq. (TjQ). Surprisingly, the long-range behavior of the PES-PS is quite different 
from ours and the one of PES-SM, and probably inaccurate due to the matching procedure 
between the long and short distances used by these authors. 



IV. VIBRATIONAL STATES OF H 3 AND D 3 

As mentioned in the introduction, one of the motivations of this study is to describe 
low-energy collisions between H 2 and H~, and to investigate the formation of H3 in such 
collisions. Therefore we need an approach which treat both bound and continuum states 
of the triatomic anion (including rearrangement of nuclei). We use the Smith- Whitten 
hyperspherical coordinates [54] : a hyper-radius p and two hyperangles 6, <f), which can be 
defined for the three identical particles by the three internuclear distances r i; (i = 1, 2, 3) as 

r< = 3- 1 /4 p ^ 1 + sin ^ sin (0 + e .) ) (5) 

where e\ = 2tt/3, e 2 = — 27r/3, and e 3 = 0. We also employ the adiabatic separation 
between hyperradius and hyperangles, which are known to be well adapted to atom-molecule 
inelastic and reactive scattering involving identical particles. The d yna mics is treated within 



the framework of the slow variable discretization (SVD) method [55145 7| that allows us to 
account easily for non-adiabatic couplings between hyperspherical adiabatic channels. We 
briefly recall below the main steps of our approach, which is discussed in greater details in 
Refs. 
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The eigenstates \I/ of three particles interacting through a potential V depending only 
on the three internuclear distances are obtained by solving the Schrodinger equation with 



the fol 
frame 



owing Hamiltonian expressed in hyperspherical coordinates p, u in the center-of-mass 



l_ _ 5 d_ 5 d_ A 2 
dp^ dp 2pp 2 



+ V 



(6) 



where p = m/y/3 is the three-body reduced mass and m is the mass of each of the three 



identical particles. The operator A above is the grand angular momentum 58 



It 



depends only on the set u of five angles, which include the three Euler angles (for the 
orientation of the molecular frame in the lab frame) and the hyperangles 9 and 0. If the 
total angular momentum J = 0, A depends only on the two hyperangles. The explicit form 
of A 2 is given, for example, by Eq. (27) of Ref . 59] : 

4 d . d 4 d 2 



+ 



24 



d . d 4 



(7) 



+ 



24 



+ 



T 2 



+ 



cos 6Jy d 



1 - sin 9 ' 1 + sin 9 ' sin 2 9 sin 2 9 d(f) ' 
where Jx, Jy, an d Jz are the components of the angular momentum along the principal 
axes of inertia. The orientation of the axes is approximately indicated in Fig. [TJ After 
rescaling the wave function \l/ as ^ = p~ 5//2 ^, the Hamiltonian for the new function tp is 
written as 



H = T p + if ad , where 



T, 



i o 2 



H, 



ad 



2p dp 2 
A 2 + 15/4 



and 



2^p 2 



+ V. 



(8) 
(9) 

(10) 



The eigenfunctions ip(p, oj) of the above Hamiltonian are sought as an expansion over the 
basis functions y a j(p,u) with unknown coefficients c aj -: 

$(j>,u) = ^2ya,j(P,u)c a ,j (11) 

a,j 

The basis functions y a j(p,u) are constructed as products 

yaj(p,w)=irj{p)<Paj(w) (12) 
where the functions 7Tj(p) are DVR-like functions localized at DVR grid points pj along the 



hyper-radius. As in our earlier study 



here we used the plane wave DVR functions. 
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The functions ip a j (u) are the adiabatic hyperspherical states obtained by solving the three- 
body Schrodinger equation at fixed p = pj, i.e. they are eigenstates of at fixed p = pj 
with eigenvalues U a (pj) 

H p = p ^ a ,{uj) = U a {p 3 )^M. (13) 

The functions U a (pj) are usually referred to as hyperspherical adiabatic potentials. In- 
serting the expansion of Eq. (TTTT) into the Schrodinger equation Hip = Eip reduces the 
equation to a generalized eigenvalue problem for coefficients c aj - with eigenvalues E 



E 



Cj'a' 



\ ~ 2p'dfP^ i ) ° a ' j '' aj + MU'WW/pta'a 

E^2(^j'\^j)pOa'j', a jCfa' , (14) 



where O a y <a j are the overlap integrals (in the u space) between adiabatic states yv,j' and 

Oa'j',aj = (Va',j'\^a,j)iu ■ (15) 

The subscripts p and u at kets in the above expressions refer to the integration coordinate 
of the bracket. 

The system (TT4"]) of equations resembles to the system of coupled-channel equations, where 
non-adiabatic couplings (w| J^|y?a)^ and ((p a >\jfw\(p a } are replaced with the overlap matrix 
elements O a ui a j. The use of overlap matrix elements instead of the derivatives of adiabatic 



56 



states with respect to p simplifies significantly the numerical solution of the equation [f 
and is the main advantage of the SVD method. 

We restricted the present computations to energies and wave functions of bound states 
for two and Dj, with total angular momentum J = 0. The hyperangle 9 varies in 
the interval [0,|] while the full interval of variation for the second hyperangle, 0, is [0,27r). 
However, for three identical particles, the interval along can be restricted to [— | ,— | ] for 
wave functions of the A[ or A' 2 irreducible representations (irreps in the following) of the 
molecular symmetry group (D 3 h), and to the interval [— f ,+f] for wave functions of the E' 
irrep. Irreps with odd parities, A'(, A%, and E" are not allowed for J = 0. This is because 
the inversion applied to the rotational part of the total wave function is reduced to the 



rotation of the system about the Y axis by 7r 
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62] . The J = rotational states are isotropic, 



so only the even parity is allowed. The hyperradius p can vary in the interval [0,oo). In this 
study, it varies from 1 a.u. to 120 a.u. for Hj a.u. and from 1 a.u. to 80 a.u. for Dj. In the 
numerical calculation, we have used equal masses m = 1837.3621 a.u. for all three atoms 
in the Hg - molecule. This value is the sum of the hydrogen mass and one third of electron 
mass. Similarly, we took m = 3670.8162 a.u. (which is the sum of the deuterium mass and 
one third of electron mass) for Dg". 

The resulting adiabatic potentials of A\ and A' 2 irreps are shown in Fig. [6] for rig and in 
Fig. [7] for Dg . The curves of E' irrep are not shown in Figs. [6]and[7J The lowest E' curves 
are almost identical to the A[ and A' 2 curves and would be indistinguishable from them in 
the figures because for low vibrational levels, the energies of levels are independent on the 
symmetry (A^A^, or E') of wave functions with respect to the proton exchange between 
the dimer and the H~ ion. (We assume here that the dimer is in a particular rovibrational 
state such that the proton exchange is made without changing the dimer state. The energy 
does depend on the H 2 rovibrational state). At large p, all adiabatic curves dissociate into 
an atom+dimer system characterized by the rovibrational state (t>d, j) of the dimer. 

Some of the adiabatic curves in Figs. M and [7J exhibit avoided crossings at energies 
above ~ 4000 cm -1 : Dynamic coupling between hyperspherical adiabatic states is mostly 
determined by the avoided crossings and responsible for nuclei exchange above the potential 
barrier identified in Fig. |2j Such transitions are much less probable at smaller energies as 
each adiabatic state is only weakly coupled to other adiabatic states of the same irrep. As 
a consequence, only one component (p a in the expansion of Eq. (fTTj) is dominant. Each 
adiabatic state for a given irrep T is correlated with a definite pair of quantum numbers 
(vd, j), and can be approximately characterized by these two quantum numbers. The value 
Q of the projection of the H2 angular momentum j on the axis connecting the dimer with 
the atom is also a relevant quantum number. For J = 0, Q is always zero and, therefore is 
not specified in Figs. [6]and[71 

To characterize completely a bound state of the trimer, an additional quantum number v t 
is needed for the excitation within each adiabatic state (along the hyperradius). Therefore, 
the bound states are characterized by four approximate quantum numbers Q, j, Vt, Vd and two 
exact quantum numbers J and V. At high energies, the mixing between different Q,j,Vt,Vd 
for given J and V becomes important. 

We summarized in Table [III the present Hg bound state energies and those of Ref. [h]] 
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FIG. 6: (Color online) The lowest hyper spherical adiabatic curves (U a (p) in Eq. (fl3l) ) for J = 
of A\ irrep (solid black lines) and of A 2 irrep (dashed red lines). The curves of the E' irrep are not 
displayed, as they are very similar to the A[ and A 2 curves for energies smaller than 4000 cnr 1 
above the lowest dissociation limit H2(0,0) + H~, and would be indistinguishable on this scale. A 
couple of avoided crossings located around 4000 cm -1 are marked. Each curve is labelled at large 
R with the pair of indexes (vd,j) corresponding to the rovibrational state of the H2 molecule. The 
energy origin is set to the lowest dissociation limit H2(0, 0) + H~. The closed circles on the (0,6) 
curve exemplifies the density of grid points in R used in the calculations. 



labelled with the set of quantum numbers JJ,£l } v t} Vd,T. In agreement with Ref. [10|], we 
found four vibrational energies levels of the A[ irrep and one other level of the A' 2 irrep. 
We also found four more vibrational levels of the A' 2 irrep. Although their energies are 
located above the Vd = 0,j = dissociation threshold, they are stable because they cannot 
dissociate due to the symmetry restriction. Their binding energies are given in the table 
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FIG. 7: (Color online) Same as Figure [6] for the D 3 molecule. The energy origin here is set to the 
lowest dissociation limit D2(0, 0) + D~. A few avoided crossings are indicated by circles. 

with respect to their first allowed dissociation limit = 0,j = 1. It is worth noticing 
that our binding energies for all levels but one are larger than the energies of Ref. 10]), 
expressing that the well depth of our potential surface is found about 17 cm^ 1 deeper than 



in Ref. 



10]). The last bound level of the A[ irrep has a very small binding energy and 



may not be fully converged. The bound state energies of D J are given in Table IIHI Since 
the Dg molecule is heavier than Hj, it has more vibrational levels. No data from previous 
calculation is available for D7. 
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Energies 


above D as 


Binding 


; energies 




Present study 


Calculations of [10] 


Present study 


Calculations of [10] 


o,o,o,o,o,^i 


2103.3 


2100.6 


-70.7 


-68.4 


0,0,0,1,0,^; 


2148.4 


2147.1 


-25.6 


-21.9 


0,0,0,2,0,^; 


2168.7 


2165.9 


-5.4 


-3.1 


0,0,0,3,0,^ 


2174.1 


2168.7 


-0.01 


-0.3 


D 00 a) 


2174.1 


2169.0 






0,0,1,0,0,^2 


2140.3 


2140.8 


-152.1 (-33.7) 


n/a (-28.2) 


0,0,1,1,0,^2 


2215.1 




-77.4 




0,0,1,2,0,^ 2 


2259.3 




-33.2 




0,0,1,3,0,^ 


2281.9 




-10.6 




0,0,1,4,0,^2 


2290.9 




-1.6 




An 6 > 


2292.5 


n/a 







a > Asymptotic energy of H +H2(u<j = 0, j = 0) dissociation 
b > Asymptotic energy of H~+H2(u<j = 0,j = 1) dissociation 



TABLE II: Comparison of energies (in cm 1 ) of H 3 bound levels obtained in the present study 
with those of Ref. |10l. Energies in the second and third columns are given with respect to Das (see 
Table H|) . The binding energies (fourth and fifth columns) are given with respect to the dissociation 
limits Dqq (for A\ levels) and Dq\ (for A' 2 levels). The energies in parentheses are given with 
respect to L>oo, to compare with Ref.jlOl]. 

V. PERMANENT DIPOLE MOMENT SURFACE FOR THE H 3 GROUND 
STATE. 

The Molpro package delivers also the components of the permanent dipole moment to- 
gether with the PES. For a fixed geometry, the obtained components of the dipole moment 
correspond to the coordinate system with the X, Y, Z axes along the principal axes of inertia 
of the molecule. In our notations, the F-axis is orthogonal to the plane of H^, therefore 
Dy = 0. For 7 = and 90° there is only one non-zero component (Dz in our notations). 
Dz becomes also equal to zero at equilateral geometries. For other geometries (7 7^0 or 
90°) there are two non-zero components, Dz and Dx- We used the following convention to 
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Energies above D as 


Binding energies 


0,0,0,0,0,Ai 


1416.5 


-126.2 


0,0,0,1,0,^ 


1474.1 


-68.6 


0,0,0,2,0,^1 


1510.1 


-32.6 


0,0,0,3,0,Ai 


1530.2 


-12.5 


0,0,0,4,0,Ai 


1539.8 


-2.9 


0,0,0,5,0,Ai 


1542.7 


-0.02 


D o a) 


1542.7 




0,0,1,0,0,A' 2 


1419.9 


-182.3 


0,0,1, i,o,a 2 


1485.0 


-117.2 


0,0,1,2,0,A 2 


1532.0 


-70.3 


0,0,1,3,0,A 2 


1563.8 


-38.5 


0,0,1,4,0,A 2 


1583.7 


-18.7 


0,0,1,5,0,A 2 


1595.1 


-7.2 


0,0,1,6,0,A 2 


1600.6 


-1.7 


0,0,1,7,0,A 2 


1602.3 


-0.01 


A)i 6) 


1602.3 





a ) Asymptotic energy of D +D2(i>d = 0,j = 0) dissociation 
b > Asymptotic energy of D~+D2(vd = 0,j = 1) dissociation 

TABLE III: Energies (in cm -1 ) of D3 bound levels obtained in the present study. Energies in the 
second column are given with respect to D as (see Table HI D as being the same for H3 and D3 ). 
The binding energies in the third column are given with respect to the dissociation limits Dqq for 
the A'i levels and Z?oi for the A' 2 levels. The -Doo and Dqi dissociation D2+D - limits are calculated 
with respect to D as and also given in the table. 

label the components: at large values of R, the component Dz is the largest of the two in 
magnitude and negative, Dx being the smallest (in magnitude) and positive. The choice of 
negative Dz at large R corresponds to the Z axis oriented from the center of mass of the 
molecule towards the H~ ion. For small values of R (for fixed 7 and r) the two components 
Dz and Dx become comparable in magnitude. For such geometries they are identified by 
ensuring a smooth variation of the components with decreasing R. Notice that because the 
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R (a.u.) 

FIG. 8: The figure shows the Dz component of the permanent electric dipole of Hg~ for one value 
of 7 = 30° and several values of r from r =0.8 (the largest value of Dz at small R) to 2.4 (the 
smaller value of Dz at small R). The ab initio values and grid points are indicated by circles. 

principal axes of inertia for Hg and Dg are the same, the obtained PDMS are the same for 
the two species. 

Figures M and [9] show ab initio values of the Dz and Dx PDM components for 7 = 30° 
for all nine values of r = 0.8, 1, • ■ ■ , 2.4 a.u. as a function of R. As one can see from the 
figures, the sampling grid of ab initio geometries is dense enough to perform an interpolation 
procedure in order to calculate Dz and Dx at any arbitrary geometry. Therefore, in the 
central region in Fig. [H we interpolate the Dz and Dx PDMS using the same 3D B-spline 
procedure as for the PES interpolation. 

In region II, when R — > 00, the analytical behavior of the largest component Dz is known: 
It decreases with increasing R as (in a.u.) Dz = er cm _n = —2R/3, where r cm _n is the distance 
between the center of mass of Hg" and the H~ ion, and e = 1 is the electron charge. The ab 
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R (a.u.) 

FIG. 9: The Dx component of the permanent dipole of for 7 = 30° and different r. The ab 
initio values and grid points are indicated by circles. 

initio Dz values confirm this behavior, so the same formula is used for extrapolation of Dz 
in region II. Outside the central region in Fig. fTJ the PDMS (except the Dz components in 
region II) are extrapolated using an empirical analytical formula. The empirical formula is 
obtained by inspecting the PDMS variations close to the boundaries of the ab initio region. 
The smallest component Dx varies at large R j)/R, where k x (r, 7) is fixed by the 

value of the PES at the final point Rf = 56.4227 a.u. of the grid, i. e. k x = V(Rf\ r, j)Rf- As 
previously, we obtain a surface for the k x (r, 7) function, which is computed at any geometry 
using a 2D B-spline interpolation. For the extrapolation of the PDMS in region I we found 
(empirically) a quadratic (for Dz) and linear (for Dx) dependencies along r with coefficients 
depending on R and 7, which are evaluated from the corresponding boundary values of Dz 
and Dx respectively. Similarly, in region III and IV, we assumed a linear dependency along 
r and R respectively, with coefficients fixed by the boundary values of the PDMS. 
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VI. SUMMARY AND CONCLUSIONS 



In the present study we have obtained the accurate potential energy surface and compo- 
nents of the permanent dipole moment for the H3 van der Waals molecule. The surfaces 
were calculated on a dense grid of geometries that covers short, intermediate, and long-range 
regions. In total, the ab initio calculations were made for 3024 geometries. The large basis 
and grid used in the calculations allows us to suggest that the obtained PES is more accurate 
than the results of the previous study jiol . [ill ] . Comparison of the long-range behavior of the 
obtained PES with the expected analytical behavior of the PES confirms this conclusion. 
No previous data on the dipole moment of is available. 

The obtained ab initio values for potential energy surface and the dipole moment compo- 
nents were used to construct Fortran interpolation/extrapolation subroutines that calculate 
the energies and dipole moments for any arbitrary geometry. The subroutines interpolate the 
surfaces using B splines inside the box of ab initio geometries and extrapolate the surfaces 
outside of the box using analytical formulas based on the theoretical asymptotic behavior. 
The subroutines are available from the journal's EPAPS service. The energy surface can 
be used for all isotopologues of , the dipole moment surfaces in the present form can 
only be used for Hg - and Dj isotopologues. For the H2D - and D2H - molecules, the dipole 
moments should be transformed to account for the different orientation of axes of inertia. 
Using the new potential surfaces, we have calculated the bound states for the Hj and Dj 
isotopologues. 

A relatively large magnitude \D\ ~ 4 a.u. of the dipole moment near equilibrium positions 
for bound vibrational states and a large size of the electronic clouds of H~ suggest that the 
cross-section for the formation of Hj stable molecules by the radiative association between 
H 2 and H~ is significant. A relatively large dipole moment and existence of several bound 
levels suggests also that the H3 can be detected using the IR photoabsorption spectroscopy. 

The calculated energies of bound states can also be used to search for in the 
cold interstellar medium with a large fraction of ionized hydrogen (to have enough free 
electrons). We have developed a model for the formation of in radiative association 
collisions between H2 and H~. In the model, H~ is formed by dissociative attachment of the 
electron to H2. If a photoabsorption signal from (in mm range) is detected, this would 
also be a signal for the presence of H~ in the ISM: H~ itself cannot be detected directly. 
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The details of the model as well as the calculated rates of radiative association collisions 
between H 2 and H~ will be discussed in a forthcoming publication. 
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